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Molten Salt Reactor (MSR) was selected as one of the six innovative nuclear reactors by the Generation IV 
International Forum (GIF). The circulating-fuel in the can-type molten salt fast reactor makes the neutronics 
and thermo-hydraulics of the reactor strongly coupled and different from that of traditional solid-fuel reactors. 
In the present paper, a new coupling model is presented that physically describes the inherent relations between 
the neutron flux, the delayed neutron precursor, the heat transfer and the turbulent flow. Based on the model, 
integrating nuclear data processing , CAD modeling, structured and unstructured mesh technology, data analysis 
and visualization application, a three dimension steady state simulation code system (MSR3DS) for the can-type 
molten salt fast reactor is developed and validated. In order to demonstrate the ability of the code, the three 
dimension distributions of the velocity, the neutron flux, the delayed neutron precursor and the temperature 
were obtained for the simplified MOlten Salt Advanced Reactor Transmuter (MOSART) using this code. The 
results indicate that the MSR3DS code can provide a feasible description of multi-physical coupling phenomena 
in can-type molten salt fast reactor. Furthermore, the code can well predict the flow effect of fuel salt and the 
transport effect of the turbulent diffusion. 
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I. INTRODUCTION 


The first molten salt reactor (MSR) was developed in late 
1940s as the aircraft propulsion at Oak Ridge National Labo- 
ratory (ORNL). The Aaircraft Reactor Experiment (ARE) [1] 
in 1954 operated successfully. In 1965, the Molten Salt Reac- 
tor Experiment (MSRE) [2] goes critical, and after six months 
of successful operation the enriched U-235 was removed and 
replaced by denatured U-233 as fissile fuel. In October 1966, 
the MSRE was the first reactor that reached criticality with U- 
233. A detailed 1000 MWe engineering conceptual design of 
a Molten Salt Breeder Reactor (MSBR) [3] was developed in 
1970s. Even though the concept looked promising, the stud- 
ies were stopped in 1976. MSR has the potential of meet- 
ing the goals of Generation IV reactors and high-level waste 
transmutation program, being better than solid fuel reactors. 
Consequently, many countries are interested in developing new 
concept MSRs, such as the FUJI series [4], Actinides Molten 
Salt TransmutER (AMSTER) [5], MOSART [6], Fast Spec- 
trum Molten Salt Reactor (FS-MSR) [7], European Molten 
Salt Fast Reactor (MSFR) [8], and Small Mobile Molten Salt 
Reactor (SM-MSR) [9]. In 2010, Chinese Academy of Sci- 
ences restarted the thorium molten salt reactor project and es- 
tablished the Center for Thorium Molten Salt Reactor Sys- 
tem (CTMSRS) in Shanghai Institute of Applied Physics, and 
CTMSRS finished the conceptual design of a Liquid Fuel Tho- 
rium Molten Salt Reactor (LF-TMSR). 

The Can-type Molten Salt Fast Reactor (CMSFR) can 
be employed to consume actinides from light water reac- 
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tor (LWR) fuel or, alternatively, to extend fissile resource 
availability through U/Pu and Th/U breeding. CMSFRs are 
highly flexible and can be configured into modified open or 
full-recycle configuration. Since 2005 the GIF has focused 
on CMSFR from thermal spectrum molten salt reactor. Many 
works have been conducted to investigate the complex behav- 
ior of the CMSEFR. Wang et al. [10] extended the SIMMER-III 
code for simulating the MOSART with extra thermo-hydraulic 
and neutronic modules at two-dimension (2D) axial-symmetric 
geometry. Nicolino et al. [11] developed a new approach to 
describe the strong coupling between neutronics and thermo- 
fluid dynamics with particular focus to the MOSART molten 
salt fast reactor at 2D axial-symmetric geometry. Zhang et 
al. [12] provided a theoretical model of flow coupling the heat 
transfer and neutronics models at 2D axial-symmetric geome- 
try and calculated of the steady characteristics of a molten salt 
fast reactor without graphite moderator in the core. 

In all the studies, the transport effect of turbulent diffu- 
sion wasn’t taken into account for the delayed neutron pre- 
cursor (DNP) concentration. In this paper, we present a new 
model to describe the inherent relations between the neutronics 
and the thermohydraulics. Based on this model, a three dimen- 
sion steady state simulation code (MSR3DS) for the molten 
salt fast reactor is developed and validated. Using the multi- 
physical coupling code, the three dimension distributions of 
the velocity, neutron fluxes, DNPs and the temperature are ob- 
tained for a simplified MOSART. 


Il. MATHEMATICAL MODEL 


More recently, several fluoride molten salt fast reactors have 
been proposed as part of the Gen IV program. These include 
the European MOSART reactor using LiF/NaF/BeF2/(TRU)F3 
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TABLE 1. Main parameters of the core [6] 


Parameters Values 
Thermal power(MW) 2400 
Reflector thickness(m) 0.2 
Inlet temperature(K) 873.15 
Inlet velocity(m - s~') 0.5 
Resident time out of the core(s) 3.94 


as a fuel salt, and an European MSFR concept using 
LiF/NaF/(TRU)F3 or LiF/NaF/(U+Th)F4. In this paper, in or- 
der to establish the theoretical model and develop the simu- 
lation code, a core configuration shown in Fig. | is adopted 
based on the core configuration of the MOSART. The main 
parameters of the core are list in Table 1. The mass propor- 
tions at equilibrium in finite critical core with 20 cm graphite 
reflector are given in Table 2. Table 3 lists six group delayed 
neutron fractions and the precursor decay constants. 


Fig. 1. (Color online) Geometry of the reactor core (Unit/mm). 


A. Thermohydraulics model 


The fuel salt flow in the molten salt fast reactor without 
moderator in the core like MOSART and MSFR is considered 
as turbulent flow. Their Reynolds numbers in the core region 
can be higher than 10°. Many models can be used to describe 
turbulent flow. The most accurate method is the direct numer- 
ical simulation (DNS), which computes the mean flow and all 
turbulent velocity fluctuations. But the DNS model is highly 
costly in terms of computing resources, so the method is not 
used for industrial flow computations. The large eddy simula- 
tion is an intermediate form of turbulence calculations which 
tracks the behavior of larger eddies. The effects on the re- 
solved flow due to the smallest, unresolved eddies are included 
by means of a so-called sub-grid scale model in the large eddy 
method. So, the demands on computing resources in terms of 
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TABLE 2. Mass proportion at equilibrium in critical core [6] 


Compositions Atomic number density/107“"cm7 
Li-6 3.856E-06 
Li-7 4.301E-03 
Na-23 1.774E-02 
Be-9 8.030E-03 
F-19 3.856E-02 
Np-237 5.713E-06 
Np-239 4.407E-11 
U-234 3.993E-07 
U-235 3.378E-07 
U-236 2.522E-07 
U-237 1.617E-09 
U-238 2.172E-09 
Pu-236 4.843E-12 
Pu-238 2.088E-05 
Pu-239 4.111E-05 
Pu-240 7.585E-05 
Pu-241 3.322E-05 
Pu-242 4.954E-05 
Pu-244 1.049E-11 
Am-241 3.728E-06 
Am-242m 1.266E-07 
Am-243 1.789E-05 
Cm-242 7.802E-07 
Cm-243 1.359E-07 
Cm-244 2.849E-05 
Cm-245 1.286E-05 
Cm-246 9.889E-06 
Cm-247 3.473E-06 
Cm-248 1.202E-06 
Bk-249 2.336E-07 
Cf-249 2.878E-07 
Cf-250 3.215E-07 
Cf-251 2.169E-07 
Zr-93 1.121E-06 
Nd-143 1.252E-06 
Nd-145 8.960E-07 
Nd-147 3.608E-08 
Pm-147 3.011E-07 
Sm-149 2.055E-07 
Sm-150 2.402E-07 
Sm-151 1.554E-07 
Sm-152 1.736E-07 
Eu-153 1.966E-07 
Eu-154 9.447E-08 
Eu-155 8.558E-08 
Gd-157 2.378E-08 
B-10 5.282E-07 


TABLE 3. Delayed neutron fractions and the precursor decay con- 
stants [12] 


Groups A(s 1) 110°”) 
C1 0.0124 22.3 
C2 0.0305 145.7 
C3 0.111 130.7 
C4 0.301 262.8 
CS 1.14 76.6 
C6 3.01 28.0 
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storage and volume of calculations are large. The Reynolds- 
averaged Navier-Stokes (RANS) equations are focused on the 
mean flow and the effects of turbulence on mean flow proper- 
ties. The computing resources required for reasonably accu- 
rate flow computations are modest, so this approach has been 
the mainstay of engineering flow calculations in engineering 
practice. 

In this study, the flow, heat transfer and turbulent charac- 
teristics in the core were obtained by solving the following 
three-dimensional, incompressible, steady state Navier-Stokes 
equations and turbulence model. The RANS equations with 
Boussiesq’s closure hypothesis and the standard k — £ turbu- 
lence model [13] were adopted. 

The continuity equation is: 


V-U=0. (1) 


The momentum conservation equation is: 


2 
ps(U-V)U =—V- (n + T 


+V- [(# + u) (VU + (VU)*)). 


(2) 


The transport equations for k and ¢ in the standard k — € 
turbulence model are: 


PU- Vk=V. (u+) vee — ps, (3) 
k 


and 
ppU Ve=V. GE H) Ve+ Cag G 
4 (4) 
-Ceai T> 
mee: Gens, S: VU + (VU)? 
Cz, = 1.44, Cz2 = 1.92, C, = 0.09 
op = 1.0, o: = 1.3. 


Where, U represents the mean velocity vector, and k, Y, pfs, 
L, H, and I are turbulent kinetic energy, turbulent dissipation 
rate, fuel salt density, dynamic viscosity, turbulent dynamic 
viscosity and Identity matrix, respectively. The equations con- 
tain five constants: Cz, Cz2, Cy, Ck and ce. The standard 
k — e turbulence model employs values for the constants that 
are arrived at by comprehensive data fitting for a wide range of 
turbulent flows. 

The energy conservation equation of the fuel salt expressed 
by temperature: 


PtsCp,ts V . (UTs) = = vale fs + Cp, co L) vT + Sfs. 
(5) 


The energy conservation equation of the reflector expressed 
by temperature: 


V. AT ret V Tref + Sref = 0. (6) 
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The inner heat sources of fuel salt (sts) and the reflector (Sef) 
in Eqs. (7) and (8) are calculated using the neutron fission re- 
actions. 


G 
ss = YE X (bg Ef), (7) 
g=1 
G 
Sret = (1 — 7) Ep 9“ (Gy « Esg), (8) 
g=1 


where Tj, Tref, AT,fs, AT refs Cp,fs and Pr; respectively repre- 
sent the temperature of the fuel salt and the reflector, thermal 
conductivity of the fuel salt and the reflector, specific heat ca- 
pacity of the fuel salt and turbulent Prandtl number; Fy, Øg, 
Xr g and y are energy released from each fission reaction, neu- 
tron flux for group g, fission cross-section for group g and the 
fraction of power released into fuel molten salt, respectively. 


B. Neutronincs model 


According to the basic conservation of the neutron number 
and DNP concentration in a control volume with the multi- 
group diffusion theory, the neutronics model of CMSFR can 
be derived. 

The diffusion equation of the neutron for group g [11]: 


1 A, 
ee 


= V- Dag") Veg (7, t) — Erg), t) 


G 
+ Xp,g(1 — 8) 5 (VE)r,g' (7) by (7, t) 
'=1 
I i G 
T 5 Xa,i,giCi (7, t) T 5 Xg'g (T)ġg (T, t). 
i=1 FG 


(9) 


The delayed neutron precursors are classified into six groups 
by half-life periods. The balance equation of the DNP concen- 
tration for group 2: 


OC; (7, t) F 
aE V- [UCT t)] 
_ Hi A 
=V. (p. i+ | VCi(r, J (10) 
I 
F Bi 5 (WE)rg (Tog (F, t) _ AiCi (r; t). 
g'=1 


The second items in Eqs. (9) and (10) indicate the flow effect 
of the fuel molten salt on the neutron fluxes and DNPs. The 
third item in Eq. (10) include two parts, the V - De ¿V C:(r, t) 
part represents the molecular diffusion of the DNP group 2 and 
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the V - (4 

Pfs Sa; 
the DNP group 2. In liquids and for turbulent flow, compared to 
turbulent diffusion, the molecular diffusion can be neglected. 


) VC; (7, t) part represents turbulent diffusion of 


In above equations, ug, Ci, Dn, g, Dei, Bis Bs V, Ug’+gs Xp,g> 
Xp,i,g and Sc, respectively represent neutron mean velocity for 
group g, the DNP for group 7, neutron diffusion coefficient 
for energy group g, DNP diffusion coefficient for group 2, the 
fraction of delayed neutron for group 2, the total fraction of 
delayed neutron, the average number of neutrons produced in 
energy group g, scatter cross section from energy group g’ to 
group g, the fission spectrums of prompt neutron for group g, 
the fission spectrums of delayed neutron for energy group g 
and delayed neutron group 2, and turbulent Schmidt number. 

In steady state conditions, the time-dependent item can be 
removed from Eqs. (9) and (10) and the effective multiplica- 
tion factor Ke is introduced. Therefore, the 3D steady state 
model of neutronics for molten salt fast reactor can be ob- 
tained. 


Drg l)o) 


G 
aes vVE)s g (7) bg (T) 


—V -[U¢,(7)] 


- =V : Da g(t) Voq(*) - 
g 


ok eat =) 


T 5 Xd,i,gri 
i=1 


WAI 
(11) 
UC] (Bat r | i) 
+E D Doy Ë) -AC 
~ (12) 


In addition, the estimate for the effective multiplication fac- 
tor keg, may be computed from the Eq. (13): 


koe l | | Loana 
= 


eff = , (13) 


[f[[ Eers a 


where, the superscript n and n — 1 represent the number of 
iterations. 


In this work, the steady neutronics model consisted of two- 
group (G = 2) neutron diffusion equations for fast and thermal 
neutron fluxes in the fuel molten salt and reflector, transport 
equations for six-group (I = 6) DNPs in the fuel salt. The 
neutron diffusion equations in the fuel salt: 


Ci(7) + 3 Lig! +g(7) bg" (P), 
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=V. Dn, ()V 61,68 (7) 
— ¥y186(7)O1,6(7) 


(1-8) `z (VX)¢,9’,t8(7) Pg’ ts(7) 


k 
eff JEL 


—V - [U¢s(")] 


U1 


(14) 
— V. Udiss(7)| =V. Dy,2,88(7)V o1,48(7) 


— X,145(7)b1,¢5(7) 
+ D132, 4s (Tr) O15 (r) . 


(15) 


The neutron diffusion equations in the reflector: 


V. Dra, retl? )V g, ret (7 )— Er, retl? )¢1, ret (7 = 0, (16) 


V y Da2 retl) V Q1 tet (7) a dy, 2, ret (7 )¢1, rett ) 
+ d4-+2 ref (7 gı, ref (7 ) 
= 0. 


(17) 


The transport equations for six-group DNPs in the fuel salt: 


V- [UC] =V (2s: + t) VCi;(r )| 
2 
T a D Egnos (T) ve 
— ACi(7) 
The effective multiplication factor kerf: 
koz! I i ip 2 VEr)g, s2 IV 
eft = (19) 


[ff Soest 


where, the subscript fs and ref respectively represent the fuel 
salt and the reflector; L145, Ur,2,f, Ur,1,ref and X, 2 ref repre- 
sent the removal cross section. 


C. Boundary conditions 
1. Thermohydraulics 


(1) Inlet boundary: at the inlet the velocity, the turbulent ki- 
netic energy, the turbulent dissipation rate and the temper- 
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ature are given: 


Uz = inlet, 
3 2 
k5 
— 0.75 
a T 


where, Uiniet is given as 0.5 m s7 t; the turbulence intensity 
I is 10%; C, is a constant assigned 0.09; l is length scale. 
Tinet 18 imposed at 873.15 K. 


(2) Outlet boundary: at the outlet a free outflow is assumed. A 
Neumann boundary condition is used for the velocity, the 
turbulent kinetic energy, the turbulent dissipation and the 
temperature. 


(3) Symmetry planes: the symmetry boundary is set for the 
velocity, the turbulent kinetic energy, the turbulent dissi- 
pation and the temperature. 


(4) Wall boundary: the boundary condition at the inner wall is 
treated by wall function method for the velocity, the turbu- 
lent kinetic energy, the turbulent dissipation and the tem- 
perature. In the outer wall the temperature is set as a con- 
stant (680 K). 


2. Neutronics 


(1) Inlet boundary: the fast and thermal neutron fluxes are im- 
posed to the vacuum boundaries. Due to the DNPs’s decay 
characteristics and resident time in external loop, the DNPs 
may return to the inlet. The DNPs at the inlet is: 


—ÀiT 
Ci inet = Ci ontlet ʻe , (21) 


where C; outiet is the DNP group 7 in the outlet and 7 is the 
resident time out of the core for fuel salt. 


(2) Outlet boundary: the vacuum boundary is for the neutron 
fluxes at the outlet, and the DNPs adopt the zero gradient 
boundary condition. 


(3) Symmetry plane: the neutron fluxes and the DNPs are all 
set symmetry boundary condition. 


(4) Wall boundary: in the inner wall the neutron fluxes use 
the coupling boundary condition and the DNPs are set to 
vacuum condition. In the outer wall the vacuum condition 
is for the neutron fluxes. 


D. Thermophysical properties and group constants 


Generally, the thermophysical properties in the thermohy- 
draulics and the group constant in the neutronics are both 
dependent on the temperature. Ignatiev et al. [6] gave the 
LiF/NaF/BeF2 solvent system properties including the density, 
the thermal conductivity, the viscosity and the heat capacity. 
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TABLE 4. The thermophysical properties of the graphite reflector at 
873.15 K [11] 


Properties Graphite 
Pret (Kg - m”) 1840 
Aryrer(W + m7". K7}) 31.2 
Cp,er(J kg? - KTP 1760 
prs (kg - m`?) = 2518 — 0.406T, (22) 


Ars (W -m~t Kt) = 0.0429 + 0.0009T, (23) 


u(Pa Í s) = 0.001e70:9942+1603.2/T (24) 


Cp.ts(J-kg~* - K~+) = 2090. (25) 


The thermophysical properties of the graphite reflector are 
considered as constants (Table 4). 

The MSR3DS code system uses DRAGON4 [14] code 
based on the XMAS-172 format libraries produced from the 
ENDF-VII.1 nuclear data to generate the two group macro- 
scopic cross-sections, the diffusion coefficients and the neu- 
tron group velocities under different temperatures. The rela- 
tions between group constants and the temperature are fitted 
by mean of the 5" order polynomial curve: 


Y.(T) =Ap + A(T — Twe) + A(T — Tavg)? 
+ A(T — Tag)? + Aa(T — Tag)“ (26) 
+ As(T — Tag)? 


where Y. denotes the group constants and Tavg is reference 
temperature. Aj—A; are fitting constants. 


E. Numerical method 


The finite volume method (FVM) is widely employed for 
solution of computational fluid dynamics (CFD) problems in 
engineering. The solution domain is subdivided into a finite 
number of small control volumes and the conservation equa- 
tions are applied to each control volume. And FVM can handle 
complex geometries. In this paper, the FVM is used for spatial 
discretization of all equations. 

We have discussed methods of discretizing the governing 
equations. This process results in a system of linear alge- 
braic equations which needs to be solved. For the discre- 
tised equations of the steady neutron diffusion and the heat 
transfer Gauss-Seidel iterative method is employed. The Pre- 
conditioned Bi-Conjugate Gradient Method (PBiCG) is used 
to solve the discretized equations of the DNPs, the momen- 
tum, the turbulent kinetic energy and the turbulent dissipation. 
And the Preconditioned Conjugate Gradient Method (PCG) 
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is applied in the discretized equations of the pressure correc- 
tion [15]. 

The SIMPLE [13] algorithm gives a method of calculat- 
ing pressure and velocities. The acronym SIMPLE stands for 
Semi-Implicit Method for Pressure-Linked Equations. It was 
originally put forward by Patankar and Spalding and is essen- 
tially a guess-and-correct procedure for calculation of pres- 
sure. It is an iterative method, and when other scalars are cou- 
pled to the momentum equations the calculation shall be done 
sequentially. The sequence of operations in the SIMPLE algo- 
rithm is given in Fig. 2. Fig. 3 is the program flow diagram of 
the coupled solver in the MSR3DS code. 


Yes 
STOP 


Fig. 2. The SIMPLE algorithm. 


Initialization 


| Solve the temperature equation 


Solve neutron diffusion and DNPs 
equations 


Solve momentum and turbulent equations 


Update the 
thermophysical 
properties and 
group constants 


Fig. 3. The program flow diagram in the coupled solver. 
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F. The description of the MSR3DS code 


The entire program diagram of the MSR3DS code is shown 
in Fig. 4. The MSR3DS code system consists of the pre- 
processing, coupled solver and post-processing module. The 
pre-processing module contains the CAD 3D modeling, mesh- 
ing and generation of the group constants. The CAD 3D mod- 
eling and meshing tool use respectively the SolidWorks and 
ANSYS ICEM CFD software. 


C_D 
DRAGON4 þ-- X-S 
Libs 
f Í 


E> va o aee 
Zoo y Yt 
I 


3D 


Geometry 


I Interface and common data Í 


H ANSYS 7 
SolidWorks ICEM CFD 


Fig. 4. The program diagram of the MSR3DS code. 


The DRAGON4 generates the group constants under differ- 
ent temperature for the MSR3DS code. The code DRAGON4 
is open source and contains a multi-group iterator conceived 
to control a number of different algorithms for the solution 
of the neutron transport equation. The SYBIL option solves 
the integral transport equation using the collision probabil- 
ity method for simple one-dimensional (1D) geometries (ei- 
ther plane, cylindrical or spherical) and the interface current 
method for 2D Cartesian or hexagonal assemblies. The EX- 
CELL option solves the integral transport equation using the 
collision probability method for general 2D geometries and 
for 3D assemblies. The MCCG option solves the integro- 
differential transport equation using the long characteristics 
method for general 2D and 3D geometries. 

The ParaView [16] code is as the post-processing module in 
the MSR3DS code. The ParaView is an open-source, multi- 
platform data analysis and visualization application. It can be 
used quickly to build visualizations and to analyze data using 
qualitative and quantitative techniques. The data exploration 
can be done interactively in 3D or programmatically using Par- 
aView’s batch processing capabilities. Furthermore, it is easily 
integrated into the MSR3DS code. 


Ill. RESULTS AND DISCUSSION 
A. Code validation 


The 3D TWIGL Seed/Blanket problem [17] is adopted 
to benchmark the neutronics calculation in the MSR3DS 
code. The effective multiplication factors calculated by CITA- 
TION [18] and MSR3DS are listed in Table 5, and the relative 
error is only 0.00147%. The fluxes normalized to 1 MW» are 
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TABLE 5. The effective multiplication factors of benchmark 


Benchmark kett Relative error/% 
3D TWIGL CITATION MSR3DS 
Seed/Blanket 0.885905 0.885892 0.00147 


shown in Figs. 5 and 6. The results verify validity of the model 
for neutronics presented in this study. 


—— CITATION-G1 
—— CITATION-G2 
—— MSR3DS-G1 
—— MSR3DS-G2 


X/em 


Fig. 5. (Color online) Distributions of the neutron fluxes calculated 
by CITATION and MSR3DS at (Z=79.2 cm, Y=40 cm). 


—— CITATION-G1 
13 —— CITATION-G2 
—— MSR3DS-G1 
—— MSR3DS-G2 


Neuton flux/(neutrons-cm”-s”) 


0.0 ee eT 
0 20 40 60 80 100 120 140 160 180 


Ziem 


Fig. 6. (Color online) Distributions of the neutron fluxes calculated 
by CITATION and MSR3DS at (X=40 cm, Y =40 cm). 


In order to evaluate the flow and heat transfer calculation 
in the MSR3DS code, a simple pipe case simplified from 
MOSART is used (Fig. 7). The thermophysical properties of 
the molten salt and the graphite at 873.15 K are applied in the 
process of validation. Similarly, we assume that the velocity 
in the inlet is 0.5 m/s in the validation of the flow calculation 
of the MSR3DS. The results are compared in Figs. 8 and 9. 
From the figures, the results calculated by the MSR3DS agree 
well with those calculated by the Fluent. Therefore, the flow 
calculation in the MSR3DS is valid. 

For validation of the heat transfer calculation in the 
MSR3DS, the temperature in the inlet is assumed as 873.15 K, 
the temperatures in the outlet and other outer boundaries are 


Fig. 7. (Color online) A simple pipe for validation of the flow and 
heat transfer (Unit/m). 


0.5055 
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0.4995 
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—— Fluent 
—+— MSR3DS 


Velocity/(m-s") 


Fig. 8. (Color online) Velocities along the symmetry axis. 
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0.10 
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L T 
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—— MSR3DS 


Velocity/(m-s”) 


-1.5 


-1.0 -0.5 00 05 10 15 20 
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Fig. 9. (Color online) Velocities along radial direction at the mid- 
plane. 


300 K. Figs. 10 and 11 show that the temperatures obtained by 
the MSR3DS are in accord with those obtained by the Fluent, 
indicating that the heat transfer calculation in the MSR3DS is 
acceptable in an engineering context. 
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Fig. 10. (Color online) Temperatures at the symmetry axis. 
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Fig. 11. (Color online) Temperatures at the mid-plane. 


B. Distributions of calculated physical fields in the core 


Using the multiphysical coupling code, we calculted under 
different conditions the 3D distributions of velocity, the turbu- 
lent kinematic viscosity, fast and thermal neutron fluxes, DNPs 
and the temperature in the core. 

Figures 12-13 show the 3D distributions of the velocity and 
the turbulent kinematic viscosity. Fig. 13, the turbulent kine- 
matic viscosity in the outlet is far greater than that of other 
locations. 

Figures 14-15 show the 3D distributions of the fast and ther- 
mal neutron flux without flow (a), with convective term (b), 
and with convective and turbulent diffusion term at Sc, = 0.7 
(c) in the balance equations for six-group DNPs. The turbulent 
Schmidt number is set to 0.7 as the default value in the ANSYS 
Fluent. The temperature in the core is set to 900 K under no 
flow condition. Figs. 14-15(a) and (b) show that the fuel salt 
flow has little effect on the distribution of the fast and thermal 
neutron fluxes. And the turbulent diffusion term hardly affect 
the distribution of the fast and thermal neutron fluxes compar- 
ing Figs. 14-15(b) with Figs. 14-15(c). 

Figures 16-21 display the 3D distributions of DNPs without 
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Fig. 12. (Color online) Velocity field in the core (U/m - s71). 


Fig. 13. (Color online) Turbulent kinematic viscosity in the 
core (vim? +s). 
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Fig. 14. (Color online) Fast neutron fluxs without flow (a), with 
convective term (b), with convective and turbulent diffusion term at 
Sc, = 0.7 (c) (Unit/neutron - m7?- s71). 


flow (a), with convective term (b), with convective and turbu- 
lent diffusion term at Sc, = 0.7 (c) in the balance equations 
for six-group DNPs. The convective term affects the distribu- 
tion of the DNPs significantly as shown in in Figs. 16-21(b), 
and the smaller the delay constant, the greater the influence of 
the flow. After considering the turbulent diffusion term, due 
to the transport effect of the turbulent diffusion term, the con- 
centrations of the DNPs decrease and the distribution of the 
DNPs in Figs. 16-21(c) are also changed observably, likewise, 
the smaller the delay constant, the greater the influence of the 
turbulent diffusion. 
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Fig. 15. (Color online) Thermal neutron fluxs without flow (a), with 
convective term (b), with convective and turbulent diffusion term at 
Sc, = 0.7 (c) (Unit/neutron - m~? - s~+). 
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Fig. 16. (Color online) C1 group precursor concentration without 
flow (a), with convective term (b), with convective and turbulent dif- 
fusion term at Sc, = 0.7 (c) (Unit/m~?). 
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Fig. 17. (Color online) C2 group precursor concentration without 
flow (a), with convective term (b), with convective and turbulent dif- 
fusion term at Sc, = 0.7 (c) (Unit/m~ °). 


In order to calculate the temperature distribution in the core, 
in this paper, the turbulent Prandtl number is imposed at 0.85, 
which is the same as the default value in the ANSYS Fluent. 
So, the distributions of calculated temperature with convec- 
tive term (a), with convective and turbulent diffusion term at 
Sc, = 0.7 (b) are shown in Fig. 22. The results in the fig- 
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Fig. 18. (Color online) C3 group precursor concentration without 
flow (a), with convective term (b), with convective and turbulent dif- 
fusion term at Sc; = 0.7 (c) (Unit/m~?). 
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Fig. 19. (Color online) C4 group precursor concentration without 
flow (a), with convective term (b), with convective and turbulent dif- 
fusion term at Sc, = 0.7 (c) (Unit/m~ °). 
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Fig. 20. (Color online) C5 group precursor concentration without 
flow (a), with convective term (b), with convective and turbulent dif- 
fusion term at Sc, = 0.7 (c) (Unit/m~ °). 


ure indicate that the distribution of the temperature in the core 
remains unchanged, because the DNPs affect the neutron flux 
slightly under the steady condition. In addition, because the 
transport effect of the turbulent diffusion term of the energy 
conservation equation at the outlet becomes more and more 
strong with the distribution of the turbulent kinematic viscos- 
ity, the fuel temperature at the outlet of the core decreases 
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Fig. 21. (Color online) C6 group precursor concentration without 
flow (a), with convective term (b), with convective and turbulent dif- 
fusion term at Sc, = 0.7 (c) (Unit/m~?). 


slightly in both cases. 
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Fig. 22. (Color online) Temperatures with convective term (a), with 
convective and turbulent diffusion term at Sc; = 0.7 (b). 


IV. CONCLUSION 


In order to investigate the complex behavior of the core in 
the CMSFR, in this present research, a new multi-physical 
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